Generation of photon-number entangled soliton pairs through interactions 
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Two new simple schemes for generating macroscopic (many-photon) continuous- variable entangled 
states by means of continuous interactions (rather than collisions) between solitons in optical fibers 
are proposed. First, quantum fluctuations around two time-separated single-component temporal 
solitons are considered. Almost perfect correlation between the photon-number fluctuations can 
be achieved after propagating a certain distance, with a suitable initial separation between the 
solitons. The photon-number correlation can also be achieved in a pair of vectorial solitons with 
two polarization components. In the latter case, the photon-number-entangled pulses can be easily 
separated by a polarization beam splitter. These results offer novel possibilities to produce entangled 
sources for quantum communication and computation. 

PACS numbers: 03.67.Mn, 03.67.-a, 05.45.Yv, 42.65.Tg 



Introduction Quantum-noise squeezing and correla- 
tions are two key quantum properties that can exhibit 
completely different characteristics when compared to 
the predictions of the classical theory. Almost all the pro- 
posed applications to quantum measurements and quan- 
tum information treatment utilize cither one or both of 
these properties. In particular, solitons in optical fibers 
have been known to serve as a platform for demonstrating 
macroscopic quantum properties in optical fields, such 
as quadrature squeezing, amplitude squeezing, and both 
intra-pulse and inter-pulse correlations. The develop- 
ment of quantum theories of nonlinear optical pulse prop- 
agation in the past years has opened a way to analyze the 
quantum features of fiber-optic solitons. Experimental 
progress in demonstrating various quantum properties of 
these solitons has also been reported, see Refs. Q-0 and 
references therein. 

In the field of quantum information processing and 
quantum computing, nonlocally entangled optical quan- 
tum states have been shown to be highly useful sources. 
The applications include quantum cryptography Q , tele- 
portation [3, and algorithms 0,13 ■ Following the Bohm's 
suggestion the entangled pairs were mostly realized 
in terms of discrete quantum variables, such as spin, po- 
larization, etc. However, the original gedanken experi- 
ment proposed by Einstein, Podolsky, and Rosen (EPR) 
utilized the continuous variables (the coordinate and mo- 
mentum of a particle) to argue that the quantum mechan- 
ics is incomplete 0]. I n 1992, Ou et al. used nondegen- 
erate parametric amplification to demonstrate the EPR 
paradox with continuous variables Later, Vaidman 
proposed a generalized method for the teleportation of 
continuous- variable quantum states fl2] | . Braunstein and 
Kimble analyzed the entanglement fidelity of quantum 
teleportation with continuous variables |13| . Quantum 
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teleportation of optical coherent states was experimen- 
tally realized by using the entanglement from squeezed 
states |l4j . After that, quantum- information processing 
with continuous variables has attracted a lot of interest 
as an alternative to single-photon schemes. 

In previous works, continuous-variable entangled 
beams have been generated by letting two squeezed fields 
(squeezed vacuum states or amplitude-squeezed 

fields 0]) interfere through a beam splitter, which math- 
ematically acts as the Hadamard transformation. By 
utilizing the continuous EPR-like correlations of optical 
beams, one can also realize quantum- key distributions 
and entanglement swapping jjjj. Thanks to these 
successful applications, squeezed states become essen- 
tial for generating entangled continuous-variable quan- 
tum states, and play an important role in the study of 
the quantum-information processing. 

It has been demonstrated that two independent 
squeezed pulse states can be simultaneously generated by 
using optical solitons in the Sagnac fiber loop configura- 
tion An EPR pulse source can be obtained by com- 
bining the two output pulse squeezed states by means of 
a 50:50 beam splitter. In contrast to this known method 
for achieving the entanglement, in this work we propose 
simple schemes for generating continuous- variable entan- 
gled states through soliton-soliton interactions, in single- 
mode and bimodal (two-component) systems, without 
using beam splitters. The quantum interaction of two 
time-separated solitons in the same polarization is de- 
scribed by the quantum Nonlinear Schrddinger Equation 
(NLSE) , and in the bimodal system, including two polar- 
izations, it is described by a system of coupled NLSEs. 
The photon-number correlation between the two soli- 
tons can be numerically calculated by using the back- 
propagation method |l6| . In addition to the transient 
multimodc correlations induced by cross-phase modula- 
tion we also find nearly maximum photon-number en- 
tanglement in the soliton pair. By controlling the initial 
separation of the two solitons, one can achieve a positive 
quantum correlation with the correlation parameter close 
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to 1. 

The single-mode system Neglecting loss and higher- 
order effects, which are immaterial for the experimentally 
relevant range of the propagation distance z, temporal 
solitons in optical fibers are described by the NLSE in 
the normalized form, 

iU z + X -U tt + \U\ 2 U = 0, 

where t is the retarded time [Tg . The input profile of the 
soliton pair is taken as 



U(z, t) = sech(z, t + p) + 7 sech(,z, t — p)e 



Hi 



(1) 



with 7, 9, and 2p being, respectively, the relative ampli- 
tude, phase and separation of the solitons. If 9 = (the 
in-phase pair), the two solitons will collide periodically 
in the course of the propagation. Otherwise, they attract 
each other within a short propagating distance, and then 
move apart due to repulsion between them. 

One can evaluate the multimode quantum fluctuations 
around the solitons by solving the linearized quantum 
NLSE 18]. For the case of two-soliton collisions Q, 
Konig et al. used an exact classical solution for the 
in-phase soliton pair, and had found that the colliding 
solitons carry both intra-pulse and inter-pulse photon- 
number correlations. However, the photon-number corre- 
lation between the colliding solitons is transient, i.e., the 
inter-pulse correlation vanishes after the collision. Unlike 
the case of the collision, the photon-number correlations 
caused by the continuous interaction between the soli- 
tons belonging to the pair initiated by the configuration 
||TJ persist with the propagation. 

In Fig. ^ we display the result of evaluation of the 
time-domain photon-number correlations for the out-of- 
phase (6 = tt/2) two-soliton pair. The correlation coef- 
ficients, which are defined through the normally ordered 



covariance, 



Cij = 



An, Ar 



Anf Ari| 



(2) 



were calculated by means of the above-mentioned back- 
propagation method [l6|. In Eq. Ahj is the photon- 
number fluctuation in the i-th slot Atj in the time do- 
main, 



Ahj 



dt[U(z, t)AU\z, t) + U*(z, t)AU(z, t)], 



At, 



where AU(z, t) is the perturbation of the quantum-field 
operator, U(z,t) is the classical unperturbed solution, 
and the integral is taken over the given time slot. As 
could be intuitively expected, nonzero correlation coeffi- 
cients are found solely in the diagonal region of the spec- 
tra (intra-pulse correlations) if the interaction distance 
is short, as shown in Fig. Qfa). As the interaction dis- 
tance increases, inter-pulse correlations between the two 
solitons emerge and grow, as shown in Figs, ^b) and (c). 
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FIG. 1: The pattern of time-domain photon-number corre- 
lations, dj, of two interacting out-of-phase solitons, with 
6 = tt/2, p — 3.5, and 7 = 1.0 in Eq. The propagation 

distance is z = 6 (a), 30 (b), and 50 (c), in the normalized 
units. The width of the time slots is At — 0.1. Note the 
difference in the bar-code scales in the panels (a) and (b), (c). 
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FIG. 2: The photon-number correlation parameter C12 for the 
soliton pair with different values of the separation (p = 3.0, 
3.5, 4.0, while 6 = 7r/2 and 7 = 1.0) in (a), and different values 
of the relative phase (0 — 0, rr/4, rr/2, while and p = 3.5, 
7 = 1.0) in (b). The inset in (a) shows the evolution of the 
interaction solitons by means of contour plots. 



In addition to the time-domain photon-number corre- 
lation pattern, we have also calculated a photon-number 
correlation parameter between the two interaction soli- 
tons, as 



ANiAN 2 :) 
AN 2 AN 2 



C12 



Here ANi^ are the perturbations of the photon-number 
operator of the two soliton, the solitons being numbered 
(first and second) according to their position in the time 
domain. 

In Fig. 121(a), we show the coefficient C12 for the soliton 
pair with the initial relative phase 9 = 7r/2, equal am- 
plitudes (7 = 1.0), and different values of the separation 



p. At the initial stage of the interaction, the photon- 
number fluctuations are uncorrelated between the soli- 
tons, C12 ~ 0. After passing a certain distance, the 
photon-number correlations between the two solitons 
gradually increase, and the pair may become a nearly 
maximum-positive-correlated one. The propagation dis- 
tance needed to achieve the maximum positive photon- 
number correlation depends on the initial separation of 
the two solitons: obviously, the interaction between them 
is stronger when the initial separation is smaller. 

On the other hand, one can fix the initial separation 
but vary the initial relative phase. For this case, the re- 
sults are shown in Fig. E{b). Similar to the case of the 
soliton-soliton collision [4(, the photon-number correla- 
tion coefficient oscillates with the period equal to that 
of the two-soliton breather, if the solitons are, initially, 
in-phase. Note that in the case which may be regarded 
as intermediate between the in-phase and out-of-phasc 
ones, 9 — tt/4, the correlation coefficient first becomes 
negative, and then positive. 

Unless 9 = (when the two solitons form a quasi- 
bound state in the form of a breather) , the two solitons 
belonging to the initial configuration will separate as 
a result of the propagation. Therefore, the interaction 
between them eventually vanishes, and thus the photon- 
number correlation coefficient may saturate before it has 
a chance to reach the value corresponding to the total 
positive correlation, which is clearly seen in the inset to 
Fig. Eta) . 

The bimodal system The time-division entangled soli- 
ton pair considered above can be separated with an op- 
tical switch. Since the time separation between the two 
solitons is, typically, on the order of a few picoseconds, 
a lossless ultrafast optical switch will be required for the 
actual implementation of the scheme. The experimen- 
tal difficulties can be greatly reduced if another scheme 
is used, which utilizes vectorial solitons in two polariza- 
tions. The model is based on the well-known system of 
coupled NLSEs 



.dU ld 2 U 



dz 
.dV 
dz 



2 dt 2 
ld 2 V 



2 dt 2 



■ A\U\ 2 U + B\V\ 2 U 
A\V\ 2 V + B\U\ 2 V 



(3) 
(4) 



Here U and V are the fields in orthogonal circular polar- 
izations, A and B being the self-phase- and cross-phase- 
modulation coefficients, respectively, with the relation 
A : B as 1 : 2 in the ordinary optical fibers [13 ■ We take 
the following initial configuration for the soliton pair [cf. 
Eq. ©], 

U = sech(i4-ii)+sech(i-ii), V = sech(i+ii)-sech(*-ii). 

Using the methods for the analysis of the classical 
vectorial solitons developed in Refs. 0, |2(1 0], we 
calculated the respective quantum fluctuations and the 
photon-number correlators numerically. It should be 
noted the total intensity of the vectorial solitons, defined 
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FIG. 3: The photon-number correlation coefficient of inter- 
acting vectorial solitons. The inset displays the evolution of 
the ^-component of the classical field. 



in terms of the circular polarizations, remain unchanged 
during the propagation, but the intensities of the linearly- 
polarized (x- and y-) components, E x — (U + V)/y/2 
and E y — (U — V)/iy/2, evolve periodically, as shown 



in the insert of Fig. |3| In this figure, we display the 
evolution of the photon-number correlation between the 
x- and y- components of the vectorial solitons, which are 
originally uncorrelated, and then become negatively cor- 
related. Recently, Lantz et al. have showed that vectorial 
solitons in the spatial domain can also develop an al- 
most perfect negative correlation between quantum fluc- 
tuations around an incoherently-coupled soliton pair |22| . 

Conclusion We have studied the quantum photon- 
number correlations induced by interactions between two 
solitons in the time-division and polarization-division 
pairs. In the former case, using the pair with suitable ini- 
tial separation and relative phase, one can generate pos- 
itive or negative photon-number-correlated soliton pairs. 
An ultrafast optical switch will be needed to separate 
the two entangled solitons into different channels. On the 
other hand, by using the vectorial solitons with two polar- 
ization components, pairs with negative photon-number 
correlations between the solitons can be generated. For 
this simple polarization beam splitter will be suffi- 

cient to separate the two entangled solitons into different 
channels. 

Such new photon-number-correlated soliton pairs fea- 
ture unique entanglement properties, which may offer 
new possibilities for applications to quantum communi- 
cations and computation. The applications will be con- 
sidered in detail elsewhere. 
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